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Coexistent State of Charge Density Wave and Spin Density Wave in 
One-Dimensional Quarter Filled Band Systems under Magnetic Fields 
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We theoretically study how the coexistent state of the charge density wave and the spin density 
wave in the one-dimensional quarter filled band is enhanced by magnetic fields. We found that 
when the correlation between electrons is strong the spin density wave state is suppressed under 
high magnetic fields, whereas the charge density wave state still remains. This will be observed 
in experiments such as the X-ray measurement. 

KEYWORDS: one-dimensional quarter filled band, charge density wave, spin density wave, quasi-one-dimensional 
organic conductors, Pauil paramagnetic limit, strongly correlated system 
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§1. Introduction 

It is found that the system with the one-dimensional 
quarter filled band becomes the coexistent state of the 
charge density wave (CDW) and spin density wave 
(SDW) due to the interplay between the on-site Coulomb 
interaction (U) and the inter-site. Coulomb interaction 
(V) by recent theoretical worksBootf The inter-site 
Coulomb interaction plays important role of the charge 
ordering, and the ground state is the coexistent state 
of 2fc F -SDW and 4fc F -CDW,B ) where fc F is Fermi wave 
vector, k-p = n/Aa and a is the lattice constant. Further- 
more, when the next nearest neighbor and the dimer- 
ization of the energy band are considered, it has. been 
indicated that 2fc F -SDW and 2fc F -CDW coexist.ii 

Quasi-one dimensional organic conductors such as 
(TMTSF) 2 X and (TMTTF) 2 X (X=C10 4 , PF 6 , AsF 6 , 
Re04, Br, SCN, etc.) are known as the one-dimensional 
quarter filled band and exhibit many kinds of ground 
state- for example, spin-Peierls, SDW, superconductiv- 
ityBS In (TMTSF) 2 PF 6 , the incommensurate SDW is 
occurred at T = 12 K, where thejrave vector is (0.5, 
0.24,-0.06) by NMR measurement MeP Recently, from the 
X-ray measurement, Pouget and Ravy argue the coexis- 
tence of 2fc F -SDW and 2fc F -CDW,l§ which, has been the- 



oretically explained by Kohayashi et alH'tf mentioned 
above and Mazumdar et al.tf 

On the other hand, in (TMTTF) 2 X, (X=Br and 
SCN), it is known that the ground state is the anti- 
ferromagnetic phase understood as Mott- Hubbard insla- 
tor phase due to the dimerization and the quarter fill- 
ing. It is clear that the wave vector of the SDW is 
commesumte, (0.5, 0.25,0l_from the measurements of 
13 C-NMREJ and iH-NMRO From thjejmgle depenence 
of satelite peak positions of 1 H-NMR,EJ'E3' the alignment 
of the spin moment along the conductive axis (a-axis) be- 
comes (t, 0, |, Oljpwhich corresponds to the recent calcu- 
lational resultsBEJ l n (TMTTF) 2 Br, 4fc F -CDW accom- 
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panied by 2fc F -SDW is found in X-ray measjjrments.ET 
This can be explained by Seo and FukuyamaEJ by using 
the extended Hubbard model. 

When the pressure is applied, the commensurate an- 
tiferromagnetic phase in (TMTTF) 2 Br at ambient pres- 
sure changes to-the incommensurate SDW phase such as 
(TMTSF) 2 X0 It is originated to the increasing of the 
hopping transfer integral, t. In the case of small t, as 
the exchange interaction is strongly influenced, the state 
becomes Mott antiferromagnetic state. When t becomes 
large, the system becomes the SDW phase due to the 
Peierls instability of the Fermi surface. The difference 
between (TMTTF) 2 X and (TMTSF) 2 X is whether the 
charge or spin ordering is localized or not. This differ- 
ence is attributed that U/t in (TMTTF) 2 X is larger than 
that in (TMTSF) 2 X. It is indicated that U/t ~ 5.0 (3.0) 
in (TMTTF) 2 X ((TMTSF) 2 X) since t in (TMTTF) 2 X 
((TMTSF) 2 X) are about-OJiiD.21 eV by the extended 
Huckel band calculations.E3 , li3 , llJ , t3 We consider the sys- 
tem in (TMTTF) 2 X ((TMTSF) 2 X) as strongly (non- 
strongly) correlated. 

In the one-dimensional system, when the magnetic 
field (H) is applied to c-axis, the amplitudes of the 
charge density or the spin moment along the c-axis in the 
CDW or SDW state by coupling of electrons with same 
spins ^ciiimpressed due to Pauli paramagnetic limit field 
(fip),BBB where H p ~ A(0)/ % /2 ( u B , A(0) is the am- 
plitude of the energy gap at H = and /ib = eh/2moc 
is the Bohr magneton. The energy band is splitted by 
Zeeman effect, so that the original wave vector at H = 
becomes the not good nesting vector. In the case of the 
magnetic field applied along the a- or 6-axis, the CDW 
and SDW are not broken, because the nesting vector is 
unchanged by Zeeman effect. In other words, the CDW 
and SDW are not influenced when the magnetic field is 
applied perpendicular to the easy axis. This picture is 
for the weak coupling system. 

In the strong coupling system, the -ff-dependence of 
the antiferromagnetic state by one-dimensional Ising 
model has been studied in the mean field approxima- 
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tionE3) The amplitude of the spin moment along the c- 
axis of the antiferromagnetic state under magnetic fields 
applied along the a-axis is easily obtained, which obey 



S z (HJ)/S z (0,j) = VI - (H/Hoy 



(1) 



where S z (H,j) is the amplitude of the spin moment 
along the c-axis at j site at H — and is the critical 
field at which the ordering of the antiferromagnetic state 
disappear.^ Since the spin of electrons are tilted to a- 
axis by the magnetic field, S z is smaller upon increasing 
H x , finally, S z becomes zero. 

In the case of magnetic fields applied to the c-axis, the 
antiferromagnetic state along the c-axis is kept, because 
if the spin is tilted to a-b plane, this tilted state is not 
unstable in weak fields. However, the paramagnetic state 
becomes more stable at higher critical field, H®, where 
H° z = H°B 

We try to analyze the coexistent state in the two cases 
when the electron correlation is strong or not, since the 
coexistent phase of CDW and SDW under magnetic fields 
dose not have been studied although the state of CDW or 
SDW under magnetic fields has been studied. When U jt 
and V/t are large (~ 5.0), we consider that the system 
is strongly correlated, because the charge and spin are 
localized as shown in Figs. 1 and 2. On the other hand, 
in small U/t and V/t (~ 1.5), the correlation between 
electrons is not strong, where there are small amplitudes 
of charge density and spin moment as shown in Fig. 7, 
8 and 9. We calculate to compare the strong coupling 
system with the non-strong coupling system by using of 
two sets of the values of U/t and V/t. 

In this paper, we calculate the self-consistent solutions 
at T = for the one-dimensional band model under the 
magnetic field perpendicular to (or parallel to) the a- 
axis based on the mean field approximation. We use the 
one-dimensional quarter filled extended Hubbard model, 
where the effect of the dimerization do not be considered 
to be simplified problems. 

§2. Formulation 

We treat the one-dimensional extended Hubbard 
model, 



n = ic + u + v, 



(2) 



it = tJ2(ClC i+1 ,„ + h.c.) - f^l^n^, (3) 



U = U ^2n tA n hl , 

i 

V = V E n ii(T n i+ha <, 



(4) 
(5) 



where C\ a is the creation operator of a spin electron at 
i site, rii >a is the number operator, g = 2, i = 1, • • • , iVg, 
iVs is the number of the total sites and a =\ and j. When 
the magnetic field is applied to (x or z)-axis (H x or H z ) 1 
Hj = H<jj (j — x and z), where H is the strength of 
the magnetic field and &j is Pauli spin matrix. In this 
model, the filling of electrons is 1/4. 



The interaction term, U and V are treated in mean 
field approximation as 

^ = E EWQ^Cfc*, i)c(k x - q, i) 

k m Q 

+ p\{Q)C\k x -QA)C{k x ,^} 

- jJ>(Q)p!(Q). (6) 



,v. 



\ ' e -iQa 



{p a {Q)C\k x ,a')C{k x -Q,u') 



+ pt,(Q)tf(k x ,a)C(k x -Q,a)} 



V 



IU ^ 

Q,a,a' 



V" e-^ Pa (Q)p* a ,(Q), 



(7) 



where / = U/Ns- The self-consistent equation for the 
order parameter p a {Q) is given by 



r (Q) =1^2 < CHk x ,a)C(k x -Q,a)> 



(8) 



We use the mean field, p a (Q), by the coupling between 
electrons with same spins. In order to simplify, we do not 
consider the case of the mean field, p~ a {Q) = I^2k < - 
C^(k x ,a)C(k x — Q,a) >, by the coupling of electrons 
with opposite spin. 

We limit the sum of the wave vector as Q = q, 2q, 3q 
and Aq (q = 2kp), because the wave vectors of 2kp = 
7r /2a and its higher harmonics should be considered due 
to the nesting of the Fermi surface in the one-dimensional 
quarter filled band. We can obtain the self-consistent 
solutions from eq. (8) by using eigenvectors obtained by 
diagonalizing K, +Um + Vm, which becomes 8x8 matrix. 
The electron density at j site, n(j), and the spin moment 
at j site, S Z (J), are given by 



JQja 



(9) 



Q 

The notation in this paper follows Seo and Fukuyama.B 
We can calculate the total energy, E, 



(ii) 



where the sum is limited to electron filling and e,^ is 
an eigenvalue. From the ordered state energy (Eos) and 
the normal state energy (-En), the energy gain (E g ) can 
be obtained by E g = Eos — ^n- 

The Pauli spin susceptibility, %, is /i|iV(0) (-/V(O) is 
the density of state i-cw— the Fermi energy and ^V(O) = 
Ns / (iirt sin akp)) and the energy gain from the 

normal state is given by —\H 2 . When E g = —xH 2 , the 
Pauli limit field, h p , is given by 



(12) 
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where h p = [iBgH p /2t. 
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Fig. 1. 2S Z as a function of v at H = 




Fig. 2. 5 as a function of v at H = 



§3. Results and Discussions 

3.1 Strong coupling 

First, we show the result at H = at U/t = 5.0, whose 
large value means the strongly correlated system. Figs. 
1 and 2 are S z and 6 as a function of V/t at U/t = 5.0. At 
< V < 0.392, the antiferromagnetic ordering ((T,T,4>4), 
i.e., S z (l) = S z {2) = -S z (3) = -S z (4)) is stabilized 
and there is no charge ordering. The spin ordering of 
(T)T,4,4) nas the wave vector of 2fcp. Above V/t = 0.392, 
the spin ordering becomes (T, 0,4,0) (S z (l) = — S z (3), 
S z (2) = S z (4) = 0) and the charge ordering (6,-6,6,- 
6) exist, where n(l) = n(3)=0.5+6, n(2) = n(4)=0.5-6, 
which can be seen in Figs. 1 and 2. These (t, 0,4,0) 
and (6,-6,6,-6) mean 2fcp-SDW and 4fcp-CDW, respec- 
tively. By including the inter-site Coulomb interaction, 
4&F-CDW is induced. 

When U/t and V/t are smaller than (~ 1.5), we under- 
stand that the system with small S z is the SDW transi- 
tion due to Peierls instability of the Fermi surface. The 



spin and charge orderings are localized if U/t and V/t 
become larger than (~ 4.0) since the amplitudes of S z 
and n are saturated. This state with the localized spin 
and charge orderings is Mott antiferromagnetic state due 
to the larger values of U/t and V/t. These are the same 
results as Seo and Fukuyama.B' 

Next, we show S z and n at H ^ and U/t = V/t = 5.0 
by using the ground state, (f, 0,4,0) and (6,-6,6,-6), in 
the strong coupling system. When the magnetic field 
is applied along x-axis (h x — ^gH x /2t), the antiferro- 
magnetic state is gradually suppressed up to the critical 
field (h c x — 2.4) and above h% the spin ordering becomes 
(0, 0, 0, 0), whereas the charge ordering is unchanged, as 
shown in Figs. 3 and 4, where 5 Z (1) = —S z (3), S z (2) = 
S Z (A) and n(l) = n(3)=0.5+<5, n(2) = n(4)=0.5-5. The 
/independence of the amplitude of S z is in good agree- 
ment with eq. (1) when we set H x as H x — 2th x / ^b<7, 
which is shown by solid lines in Fig. 3. 

In the case of h z = ^-QgH z /2t ^ 0, the alignment 
of the spin moment, (t, 0,4,0) when h% = is kept 
as h c z increases, but, the system becomes (4,0,4,0) at 
h° z — 2.4, as shown in Fig. 5. This h° z has the same value 
of h x . The charge ordering is not changed upon increas- 
ing h z , as shown in Fig. 6, where n(l) — n(3)=0.5+6, 
n(2) = n(4)=0.5-6. These (4,0,4,0) and (6,-6,6,-6) are 
4fcp-SDW and 4fcp-CDW, that is, above h z the system 
becomes paramagnetic state to stay to be localized due 
to larger U/t and V/t. 

The h x - and h z - dependences of the spin ordering in 
our results can be understood by the mean field solu- 
tions for strong coupled Ising model mentioned in the 
introduction^ 

The charge ordering is unchanged by the magnetic 
field. The magnetic field is not contributed to localized 
electrons made by larger V/t. 




Fig. 3. When U/t = V/t = 5.0, 2S Z as a function of h x . The 
solid lines are written by eq. (1). 
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Fig. 4. When U/t = V/t = 5.0, 8 as a function of h x 
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Fig. 5. When U/t = V/t = 5.0, 2S Z as a function of h z 




3.2 Non-Strong coupling 

We calculated the solutions as a function of V/t at 
U/t = 1.5. There were two solutions ((4,4,4,4) and 
(t, 0,4,0) and (5—5,5,-5)), which are shown in Figs 7, 
8 and 9. These are for the state with 2fcp-SDW and the 
coexistent state with 2fcp-SDW and 4fcp-CDW, respec- 
tively. It is found that these amplitudes of S z and 5 
are small. In the region of < V/t < 1.5, the energies 
with these solutions are nearly the same, namely, these 
states are degenerate. Therefore, we analyze h x - and 
/^-dependences of S z and n for two solutions. 

We calculate the case of U/t = V/t = 1.5 at H ^ 
0. On using (t,T,4,4) when H = 0, for h z this state is 
changed to (4,4,4,4) at h z = 0.0155, but, for h x , S z is 
unchanged, which are shown in Figs. 10 and 11. This 
suppression for h z comes from Pauli paramagnetic limit 
since h c z corresponds to h p ~ 0.0125 obtained by using 
cq. (12) and E g /N s ~ 0.0000175 calculated at U/t = 
V/t = 1.5 and H = 0. For h x , there is no Pauli limit, 
because the nesting vector of the SDW is not affected 
by Zeeman splitting. In Fig. 11, S z linearly increases as 
h z increases, which means that the paramagnetic state 
is stabilized by the magnetic field. 

For (4,0,4,0) and (5,-5,5,-5), when the magnetic field 
is applied along x-axis, the spin and charge ordering do 
not change, as shown in Figs. 12 and 13. However, for 
h z , both orderings of the spin and charge disappear at 
h z = 0.0155, which is corresponding to Pauli param- 
agnetic limit field, h p ~ 0.0125, as shown in Figs. 14 
and 15. This is due to the effect of Pauli limit, too. 
The coexistent state of 2fc F -SDW and 4A: F -CDW with 
(4,0,4,0) and (5,-5,5,-5) is changed to paramagnetic state 
with (4,4,4,4) an d (0,0,0,0). It is seen that the linear in- 
creasing of S z above h° z . 



U/t=L5 
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Fig. 7. By using of (t,T,U) at U/t ■ 
function of v at H = 0. 



1.5 and V/t = 0, 2S Z as a 



Fig. 6. When U/t = V/t = 5.0, <5 as a function of h z 



Coexistent State of Charge Density Wave and Spin Density Wave in One-Dimensional Quarter Filled Band Systems 



5 



U/t=1 .5 
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Fig. 8. By using of (t,0,|,0) at U/t = 1.5 and V/t = 0, 2S Z as a Fi S- n - When U/t = V/t = 1.5, 2S Z as a function of h z by using 
function of v at H = 0. of (T.T,i>J.) at h * = 0. 
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Fig. 10. When U/t = V/t = 1.5, 2S Z as a function of h x by using Fig _ 13 whcn jj/ t = y/ t = L5j g as a f unct ion of ^, by using 
of (t,T,U) at h x = 0. of (t,0,|,0) at h x = 0. 
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Fig. 14. When U/t = V/t = 1.5, 2S Z as a function of h z , by 
using of (T,0,i,0) at h z = 0. 
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Fig. 15. When U/t = V/t = 1.5, <5 as a function of h z , by using 
of (1,0,1,0) at h z = 0. 



3.3 Comparison with Experiments 

The coexistent state with CDW and SDW is realized 
due to the inter-site Coulomb interaction even if the cor- 
relation between electrons is strong or not. 

In the non-localized SDW system, when the magnetic 
field is applied parallel to the easy axis, both orderings 
of SDW and CDW disappear at the critical field of the 
Pauli paramagnetic limit. In the case of the magnetic 
field perpendicular to easy axis, both orderings of CDW 
and SDW are unchanged. In (TMTSF) 2 X, for example, 
since the easy axis is 6-axis, it is expected that when the 
magnetic field is applied to &-axis both orderings become 
to disorder at Pauli limit field. 

On the other hand, in the localized antiferromagnetic 
state such as (TMTTF) 2 X, the charge ordering is not 
suppressed in both cases of the magnetic field applied 
to parallel to and perpendicular to the easy axis. Thus, 
4fcp-CDW may be observed from the X-ray measurement 
even if the magnetic field is applied to a, b and c-axis. 
It is a means of finding whether the system is strongly 
correlated or not. 



In (DCNQI)2Ag, which are strongly correlated sys- 
tem such as (TMTTFkX_4fc F -CDW has been observed 
at zero magnetic fieldSEil 1 Even under high fields, the 
charge ordering should be appeared. 

§4. Conclusions 

We theoretically study the coexistent state of CDW 
and SDW under the magnetic field. As a result, in the 
case of the strongly coupling system, although the spin 
ordering is suppressed at high fields, the charge ordering 
still remains. When the coupling is not so large, the 
CDW and the SDW disappear at Pauli paramagnetic 
limit field. 

These features of the coexisitent state of CDW 
and SDW under magnetic fields should be observed 
in the strongly correlated system (non-strongly corre- 
lated system) such as (TMTTF) 2 A and (DCNQI) 2 Ag 
((TMTSF) 2 X). 
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